Complex Networks: effect of subtle changes in nature of randomness 
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In two different classes of network models, namely, the Watts Strogatz type and the Euclidean 
type, subtle changes have been introduced in the randomness. In the Watts Strogatz type network, 
rewiring has been done in different ways and although the qualitative results remain same, finite 
differences in the exponents are observed. In the Euclidean type networks, where at least one finite 
phase transition occurs, two models differing in a similar way have been considered. The results show 
a possible shift in one of the phase transition points but no change in the values of the exponents. 
The WS and Euclidean type models are equivalent for extreme values of the parameters; we compare 
their behaviour for intermediate values. 
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I. INTRODUCTION 

Ever since the discovery of small world effect in many 
natural and artificial networks, many models have been 
proposed to mimic their features [lj . Basic properties like 
the small world effect can be obtained from the Watts 
Strogatz (WS) model Q while a prototype of the scale- 
free model is the Barabasi Albert model [3[ ■ For networks 
in which links depend on the geographical distance or 
age separating the nodes, a few models have also been 
proposed [HQ- 

In all these models, an essential ingredient is random- 
ness. Randomness is involved in the way the nodes are 
connected in the network. For example, in the origi- 
nal WS model, links between nearest neighbours were 
rewired with a certain probability to a randomly chosen 
node. 

While randomness is the key feature, one may ask the 
question, to what extent does its nature affect the net- 
work? For example, one can have a variation of the WS 
model, in which links are added randomly keeping the 
original links intact @. Thus one can conceive of two 
kinds of WS models, the rewiring type and the addition 
type. The network properties are in general identical for 
the two types. 

Another class of models which has been studied exten- 
sively includes the evolving networks. In these networks, 
a new node prefers to get attached to older nodes char- 
acterised by some particular feature like degree, age or 
geographical proximity etc. In an Euclidean network, 
links are established between nodes with a probabil- 
ity which depends on the Euclidean distance separating 
them. These probabilities can be tuned to change the 
nature of the randomness. Changing the degree of the 
preferential attachment in an evolving network has been 
shown to lead to drastic changes its behaviour Q . How- 
ever, this variation is much more severe compared to the 
example of addition and rewiring type WS models. 

In this paper, we have considered two classes of models. 
One is the WS type with rewiring and the other class is 
Euclidean networks in which links are attached with the 
probability proportional to l~~ a where I is the Euclidean 
distance separating them. Within both classes, we now 



introduce subtle changes in the randomness (to be elab- 
orated in the appropriate sections) and see how it affects 
the (static) network properties. At one level, we discuss 
the effect induced by the change in randomness within 
each class and on a different level, we compare the WS 
and the Euclidean networks which are equivalent for ex- 
treme values of the parameters. 

In section II, we briefly describe the models and the 
network properties which have been studied. The WS 
type networks have been discussed in section III and in 
section IV, the Euclidean network results have been pre- 
sented. Comments and discussions are made in Section 
V. 



II. THE MODELS AND THE NETWORK 
PROPERTIES STUDIED 

To construct the WS type of networks, we start with 
a regular network with 2K — 4 nearest neighbours. In 
the first type of WS network, henceforth referred to as 
WSA, we rewind both the first and second neighbours 
with probability p. In the other type, only the second 
neighbours are rewired and we call it the WSB. 

WSB may be looked upon as a particular case of WS 
networks with 2K original links for each node and where 
rewiring of the 2nd, 3rd....i\th links are done only, i.e., 
the first nearest neighbour links are kept intact. Such a 
case was recently considered with rather high values of K 
and the results were found to be identical with WS [8]. 
This result is not surprising, as the number of rewired 
links for large K is considerably larger than the number 
of links which are not rewired. In case of K = 2, i.e., the 
present case, however, the effect (if any) of this kind of 
partial rewiring on the network properties is expected to 
be much more prominent as the number of links which 
are being rewired and those which are not, are both equal 
to L (the number of nodes), in model WSB. 

In the Euclidean network, nodes placed at a distance 
I along a one dimensional chain are linked with a prob- 
ability P(l) oc l~ a . In the first type of model which we 
call Euclidean A type (EA) , nodes at any I > 1 are linked 
while in the second type, EB, the first nearest neighbours 
are always linked while other nodes are linked with the 
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probability P(l) with I > 2. 
networks is given in Fig. [TJ 
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FIG. 1: (Color online) Schematic diagram for different net- 
work models. Average degree is 2K = 4 in each network. In 
the regular network both the first and second nearest neigh- 
bour are present. In WS network A both the neighbours are 
rewired while for WS network B only second neighbours are 
rewired with probability p. In Euclidean network A nodes at 
an y I > 1 are linked while in Euclidean network B, the first 
nearest neighbours are always linked while other nodes are 
linked with the probability l~ a with / > 2. 

In both WSB and EB, a broken network will never ap- 
pear. The Euclidean network of type B was first studied 
in the context of navigation 4] and later in relation to 
polymers [91-41 II] where short cuts or bridges appear in a 
comparable manner. Other studies on the EB were later 
conducted where issues like clustering properties, navi- 
gation with local information and phase transition of the 
Ising model were addressed (l2 - Ha |. The results of the 
various earlier studies indicate that the EB model be- 
haves as a regular network for a > 2 and like a small 
world for a < 1. For 1 < a < 2, there is some contro- 
versy regarding the behaviour of the model; while some 
studies suggest that it is small world like, there are oth- 
ers which indicate a regular latti ce ( of dimension greater 
than one)-like behaviour [l6| . In pjj , it was claimed that 
there are two transition points and the region < a < 1 
is a random graph with vanishing clustering tendency. In 
these studies, EB has been considered with two nearest 
neighbours plus one (on an average) long range link for 
each node. The main motivation for the present study 
comes from the fact that EA has never been studied ear- 
lier. Here we consider for EA, four random long range 
links for each node on an average and for EB, two near- 
est neighbour links plus two random long range links (on 
an average). This is done to make the Euclidean net- 
works comparable with the Watts Strogatz like network 
for which the average degree is four. 

The network properties we have studied are the ba- 
sic ones like average shortest paths, clustering coefficient 



and degree distribution. The average shortest paths have 
been calculated exactly using a burning algorithm. 

The two types of models, WS and Euclidean, are equiv- 
alent and identical to a random network in the limit p = 1 
and a — 0. They are also equivalent for p = and 
a = oo, corresponding to a regular chain of nodes with 
four nearest neighbour links for each node. For this reg- 
ular network the average shortest path is 
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where L is the number of nodes or system size and X is 
the remainder when (L — 1) is divided by 4. 

Clustering coefficient of the ith node is defined by 
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where Tj is the number of closed triangles attached with 
the node i and fc, its degree. The average clustering co- 
efficient of the regular network at p = or a = oo is 
0.5. 

Degree distribution is another important property of a 
network. Both the WS type and Euclidean networks are 
homogeneous, which means that the degree distribution 
will have a typical scale. For the regular network ob- 
tained in the the extreme limits a — >■ oo and p = 0, the 
degree distribution is simply a delta function at k — 4. 



III. WATTS STROGATZ TYPE NETWORKS 

We have calculated the average shortest paths as a 
function of the number of nodes (L) and rewiring proba- 
bility (p). For both WSA and WSB, the average shortest 
path varies as ln(L) for p ^ 0. But when the variation 
against the rewiring probability p is studied, there ap- 
pears to be some quantitative difference in the results 
between WSA and WSB. 

It is now widely accepted that the average shortest 
path obeys the general scaling form [1] 
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where N is the total number of nodes in a space of di- 
mension d and 

f(KpN) = for KpN»l. 

KpJSI 

In the present case, the dimension of the system d = 1 
and N = L, such that equation ([2} becomes 



(s) ~ (1/i^) 



2 Mp)+\n(KL) 



P 



for KpL » 1. (3) 



However, from the numerical results for finite systems 
shown in Fig. [2j we find the following scaling form to be 
more appropriate 
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is another important quantity in a network. We find 
that in the WS networks, the maximum degree is about 
10, the clustering coefficients show an exponential decay 
with k within this small range of k (Fig. [4j. For WSA 
the exponential decay becomes much weaker as p — > 1, 
while for WSB, the slope in a log-linear plot shows a 
much stronger dependence of C(k) on k even for large 
value of p. The clustering coefficients rapidly decrease as 
p increases, as is expected, and at p = 1 becomes very 
close to zero showing no dependence on k. The clustering 
coefficients are obviously larger in WSB than in WSA in 
magnitude (Fig. @| and do not show any dependence on 
the system size for either model. 



FIG. 2: (Color online) Variation of average shortest path as 
a function of rewiring probability p for the WSA and WSB 
models (L — 5000). Inset shows the variation of b with system 
size. 



where b is not exactly equal to unity and depends on 
the system size. The value of 6, when plotted against L, 
indicates that it saturates at large values of L, with the 
saturation values close to but not exactly equal to unity 
( see inset of Fig. [5]). Moreover, the saturation values of 
b (~ 0.95 for WSA and - 0.91 for WSB) have a finite 
difference for the two models at least over the range of 
values of L considered here. 

The average clustering coefficient for p — is 0.5 as 
mentioned before. For p > 0, two neighbors of a node 
i that were connected at p — are still neighbors of i 
and connected by an edge with probability (1 — p). Since 
there are three edges that need to remain intact for a 
cluster, C{p) ~ C(0)(1 -p) 3 in WSA. For WSB C(p) ~ 
C(0)(1— p) as we do not rewire the first neighbours in this 
case. This is exactly what we observe in the numerical 
studies (Fig. |3]). Hence WSA and WSB have clearly 
different exponents as p — > 1, the point at which the 
transition to random network is obtained. 
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FIG. 3: (Color online) Average clustering coefficient C(p) as 
a function of rewiring probability p for the WSA and WSB 
model for different system sizes are shown. C(p) is fitted to 
the form C(p) ~ (1 - p) x with x = 3 (WSA) and 1 (WSB). 

The clustering coefficient as a function of the degree 
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FIG. 4: (Color online) Average clustering co-efficient as a 
function of the degrees (C(k) vs k) for different rewiring 
probabilities for the WSA and WSB model for system size 
L = 2000. 

The degree distribution, which at p = is a delta func- 
tion, shows a spread as p is made nonzero. The variance 
of the degree distribution as a function of p can be stud- 
ied. The results show that the nature of variance is iden- 
tical in WSA and WSB, it shows a monotonic increase 
with p as expected [Fig. [5] . 

One obvious difference between the WSA and WSB 
regarding the degree distribution is that P(k < 1) = 
for WSA and P(k < 2) = for WSB because of the way 
the rewiring is done [l|. On the other hand the other 
probabilities, e.g., P(k > 2) for WSB are nonzero and 
almost independent of the system sizes considered. 



IV. EUCLIDEAN NETWORKS 

In the Euclidean network, we have similarly two mod- 
els which differ in a subtle manner as described in section 
II. 

In Fig. |6l we present the result for the shortest paths 
against L for both EA and EB models for different a. The 
log-linear plot shows that there could be a deviation from 
the behaviour (s) oc logL even below a = 2.0. This is in 
consistency with some earlier results [ll|, EH [3 which 
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FIG. 5: (Color online) Variance of the degree distribution as 
a function of rewiring probability p for the WSA and WSB 
model for system size L — 2000 are shown. Inset shows the 
degree distribution P(k) against k. 



claim that there is another phase transition at a = 1. 
We have not attempted a very detailed study to check 
exactly for which value of a the small world feature (i.e., 
(s) oc log L) vanishes as one needs huge system sizes for 
that and exact numerical calculations becomes difficult. 
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FIG. 6: (Color online) Log-linear plot of variation of average 
shortest path against system size for different values of a for 
EA and EB. The behaviour of (s) apparently deviates from 
(s) oc \ogL even below a = 2.0. 
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FIG. 7: (Color online) Comparison of average shortest path 
against a for EA and EB. In the region 2 < a < 10 we have 
finite differences in (s). 




FIG. 8: (Color online) Variation of average clustering coeffi- 
cients with q for L — 1000 and L — 2000 are shown for EA 
and EB models. Insets show the scaling plots for EA (top) 
and EB (bottom). 



pected value 1/2. 




When we plot the shortest paths against a, it is appar- 
ent that there is a transition at a — 2 beyond which both 
EA and EB behave as regular networks. This is shown 
in Fig. [7] This result is well established in all previous 
studies. When we compare the shortest paths for EA and 
EB, we find that finite differences in the results appear 
for 2 < a < 10 and the transition at a = 2 looks much 
sharper in EA. At large values of a, for both EA and EB, 
the shortest paths saturate at a value corresponding to 
that obtained from ([1} as expected. 

The average clustering coefficients calculated from 
these networks are shown as functions of a in Fig[SJ Once 
again, the values saturate at large values of a to the ex- 



6 5 

FIG. 9: (Color online) A toy model of EB type in which the 
average degree is four. The clustering coefficient of this model 
is ~ 0.53, larger than 0.5, the clustering coefficient of a regular 
network. 

There are, however, several interesting facts to be 
noted. First, we find that clustering coefficient C(a) is 
marginally larger than 1/2 for some values of a > 2 in- 
dicating that a network with some randomness can have 
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FIG. 10: (Color online) Variation of clustering coefficients 
against degree for EA and EB for L = 1000. 

larger clustering coefficients compared to the correspond- 
ing regular one. We present in Fig. |H] a toy model which 
has larger clustering coefficient compared to the regular 
one to show that it can be possible. 

In the Euclidean models, the system size dependence 
of the clustering coefficient shows an interesting result. 
Normally, one does not expect that clustering coefficient 
will have size dependence but we do get such a result here 
in these type of networks. According to PJj], it is indeed 
possible if there is a transition point below which the 
clustering is zero and above which it is finite. Hence, as 
in pjj , we attempt a data collapse by scaling the variables 
assuming C[a) to satisfy the following behaviour 

C(a) xL-"<f[(a-a c )L 5 }, (5) 

where the scaling function is a constant for very small 
values of the argument and varies as [(a — a c )i 5 ] 7 /' 5 in 
the opposite limit. 

Using the above form, we find excellent data collapse 
for different sizes with 7 = S ~ 0.2 for both EA and 
EB ( see insets of Fig. However, the value of a c 

is different; a c ~ 1.2 for EA and a c ~ 1 for EB. The 
above analysis gives us a lot of important information: 
the exponents are independent of EA and EB while the 
transition points are not. The scaling functions are also 
different for the two. The values of 7 and 5 being almost 
same showing that the clustering is a linear function of 
a — a c at a >> a c . 

As a function of the degree fc, the clustering coefficients 
are plotted against k in Fig. [TU] The values of k are 
limited and C{k) shows weak dependence on k and the 
variations do not show any drastic change with a in either 
EA or EB models. 

Next we consider the degree distributions. For both 
EA and EB, these are shown in Fig. [TTJ For EA, mini- 
mum value of the degree can be zero while in EB, there 
is a lower cutoff at 2. Consequently, the degree distri- 
butions, both of which are peaked around the average 
value 4, have different variances as function of a which 
are plotted in the inset. Here we find that the variance 



f \ 






L=2000 EA — ^ 
a=1.0 EB — x— 


1/ \ *i 


10 
















u 

§1 
















(0 










0.1 










V 01 


V 10 



k 

FIG. 11: (Color online) Comparison of degree distribution of 
EA and EB for system size L = 2000. 

for the two models are significantly different upto a ~ 2; 
for EB, it is constant upto a = 1 and falls off smoothly, 
but for EA, it has a different constant value till a = 1, 
develops a kink between a — 1 and 2 before converging 
with the EB values beyond a = 2. 



V. SUMMARY AND DISCUSSIONS 

We considered network models with a fixed average 
degree 2K = 4 where the randomness in linking the nodes 
has been conceived in different ways. 

The WSA and WSB networks, as demonstrated in sec- 
tion III, do not show qualitative differences in the main 
features. Finite differences in some exponents are there; 
e.g., the value of b appearing in the shortest path expres- 
sion shows a difference. It is expected that for WSA the 
value of b should approach unity in the thermodynamic 
limit. Similarly, for WSB too the value may approach 
unity for infinite systems. However for the finite sizes 
considered here, we get different values for WSA and 
WSB neither of which is equal to unity. Since we evaluate 
the shortest paths exactly (numerically) it is difficult to 
study very large sizes. Hence the results may have finite 
size effects, and if b really approaches unity for both WSA 
and WSB in the thermodynamic limit, one can conclude 
that the finite size effects are more for WSB. 

The average clustering coefficient vanishes as p — > 1 
in both WS types of models with an exponent that is 
markedly different and can be explained why. The vari- 
ance of the degree distribution on the other hand shows 
difference in the magnitude only. The clustering coeffi- 
cients as functions of degree are qualitatively similar in 
WSA and WSB models although in the latter, the de- 
pendence on k is stronger as p — > 1. 

For the EA and EB models also, we find the qualitative 
results to be similar. The presence of the phase transition 
at a — 2 is evident in both models, the shortest paths in 
EA showing a sharper transition. Once again, we have 
estimated the shortest paths exactly which restricts the 
system sizes considerably such that the position of the 
second transition (below a = 2), if any, is difficult to 
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establish numerically. However, analysis of the clustering 
coefficients indicate the presence of another transition at 
a c , the exact value of a c being different in EA and EB. 
The exponents, however, are identical and comparable 
to the model with a lower average degree studied earlier 

- 12 J- 

One interesting point to note is that apart from the 
shortest paths, no other quantity shows any difference 
for EA and EB beyond a — 2. Although the clustering 
and degree distribution for EA and EB are identical for 
a > 2.0, there is an interesting behaviour of the clus- 
tering coefficient which increases beyond the value 0.5 
corresponding to the regular network. 

The Euclidean model B has been considered to be ap- 
propriate to model linear polymers or self avoiding walks 
[ToL [TTj . The value of a which will be comparable to a 
linear polymer is less than 2. However, the absence of 
magnetic phase transition on a polymer chain indicates 
that it behaves as a regular network while for the Eu- 
clidean model B, the behaviour at a < 2 is definitely 
different from that of a linear chain. This indicates that 
the EB does not really mimic the behaviour of a linear 
polymer possibly because in the former, the long ranged 
links are uncorrelated which is not true for a polymer 
chain where the bridges are correlated. 

As mentioned in section II, the WS and Euclidean 
models are equivalent at extreme limiting values of the 
parameter a and p. So it also makes sense to compare 
the two when the parameters assume intermediate values. 
The main difference is of course that there is at least one 
non-trivial phase transition in the Euclidean models. In- 
creasing a beyond 2, we have an extensive region where 
regular network-like behaviour can be observed. No such 
study is possible for the WS type networks, which be- 
haves like a small world network for any nonzero value of 
p / I . In case one concludes that the Euclidean networks 
have a random nature below a c , then there is a finite re- 



gion < a < a c with random network behaviour. For 
the WS type networks, random network behaviour exists 
for p = I only. 

In both models, the degree distributions have exponen- 
tial tails. The variance in the Euclidean network however 
shows a constant value over a considerable region, which 
is not observed in the WS models. This may be due 
to the random network like behaviour of the Euclidean 
models in the entire region a < a c . 

The clustering coefficient as function of the parameter 
p or a shows that it has a scaling behaviour as the point 
p = 1 or a = a c is approached. While for WS type 
networks, the exponents are appreciably different for the 
different schemes of rewiring, for the Euclidean networks, 
the two exponents 7 and S are identical for EA and EB 
while the value of a c shows a difference. 

For WSA each node has at least K/2 = 2 edges after 
the rewiring process So there are no isolated nodes 
and the network is usually connected. But for EA there is 
no such restriction and one gets broken networks for EA 
with much higher probability. In our simulations broken 
WSA networks were not obtained at all (unless L is very 
small) while for the EA, such networks, even for large 
values of L were generated and were not considered for 
the various calculations. 

In summary, we have made a detailed study of 
different properties of some models of networks in which 
the randomness has been incorporated in different ways. 
We conclude that small changes in the randomness 
in network models do not alter the gross qualitative 
features. 
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